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Abstract: The phylogenetic relationships of four species, Bufo melanostictus, Hyla chinensis, Rana limnocharis 
and Rana guentheri, which belong to three differrent families of Anura, were detected with RAPD technique. The ge- 
nomic DNA of each species was amplifted with 19 random primers. 16 primers given clear amplified bands were used for 
analysis and the genetic distances between four species were calculated. The results show that RAPD bands obtained by 
all 16 primers evinced different degree polymorphisms. The genetic distance between R. limnocharis and R . guentheri is 
the nearest, that between B . melanostictus and H . chinensis is the second nearest, and that between B. melanostictus and 
R. guentheri is the furthest. The different distances also indicated that the relationship between Bufoidae and Hylidae is 
closer than that between Bufonidae and Ranidae at genomic DNA level. In concordant with the conclusion of the morphol- 
ogy. chromosomal and mitochondrial DNA studies. our resulta provide a new evidence of the systematic evolution of the 


three families of Anura at DNA molecular level. 
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As a method for DNA fingerprinting based on the 
polymerase chain reaction (PCR). the Random Am- 
plified Polymorphic DNA (RAPD), which developed 
by Williams et af. (1990), Welsh & McClelland 
(1990) and Welsh et af. (1991), is currently receiv- 
ing particular attention. Because of its extreme sim- 
plicity, it has been extensively used as a powerful tool 
for such tasks as gene mapping, population and pedi- 
gree analysis, phylogenetic study and bacterial strain 
identification. The basic strategy involves the PCR 
amplification of random fragments of genomic DNA with 
single or multiple primers of arbitrary sequence. Poly- 
morphism between individuals (or strains) is detected 
as differences between the pattem of DNA fragments 
amplified from the different DNAs using a given primer 
{s}. 

Anura is the biggest order in Amphibia. There 
are about 3 500 existing species in this order that dis- 
tribute widely all over the world. Many studies have 
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been working on the population genetic and systematic 
evolution of Anura by morphology, chromosomal and 
mitochondrial DNA. 

In the classification system that revised by Noble 
in 1931 (Liu & Hu, 1961), Bufoidae and Hylidae 
were in the same suborder (Procoela, Noble) while 
Ranidae in other suborder (Displaasiocoela, Noble), 
for the vertebra of Bufoidae and Hylidae are both pro- 
coelous while the front 7 vertebra of Ranidae are pro- 
coelous but the 8” vertebra is dicoelous. According to 
the probable relationship of major groups of Anura that 
suggested by Kluge & Farris (1969) based on morphol- 
ogy, the relationship between Hylidae and Bufoidae is 
closer than that between Bufonidae and Ranidae. In 
the classification system of existing Anura that suggest- 
ed by Webb et al. (1981), Bufoidae and Hylidae 
were classified into the same superfamily ( Bufonoidea ) 
while Ranidae was classified into other superfamily 
( Ranoidea) . 
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In the field of chromosome. Miura (1995) con- 
firmed by the late replication banding patterns tech- 
nique that there are homogenous chromosome among 
Hylidae, Bufoidae and Ranidae. And some regions or 
most regions in the late replication handing pallems of 
six macrochromosome between two species which re- 
spectively belong to Bufoidae and Hylidae showed ho- 
mology while only four macrochromosome between Hyl- 
idae and Ranidae showed homology. Correspond with 
the result concluded by morphologic ( Duellman. 
1975), this result manifested that in cytogenetics the 
relationship between Hylidae and Bufoidae is closer 
than thal between Bufonidae and Ranidae. 

In the field of mitochondrial DNA, the study by 
Hay et af. (1995) has analyzed most families in Am- 
phibia, included those in Anura. The result also man- 
ifested that the relationship between Hylidae and Bu- 
foidae is closer than thal between Bufonidae and 
Ranidae. These studies all manifested that the rela- 
tionship belween Bufoidae and Hylidae is closer than 
that between Bufonidae and Ranidae . 

In order to provided new evidence of the syslemal- 
ic evolution of the three families of Anura at genome 
level, we checked the phylogenetic relationship of four 
species; Bufo melanostictus, Hyla chinensis, Rana 
limnocharis and Rana guentheri which belong to the 
three differrent family of Anura, using RAPD tech- 


nique. 
1 Materials and Methods 


1.1 Experimental animals 

Referred to the report of Wang et al. (1996), 
the intra-species genetic distance is greatly closer than 
the inler-species genetic distance, so inter-species 
comparison can be made on a random individual of dif- 
ferent species. In this study, the experimental ani- 
mals are 4 random individuals respectively from 4 fami- 
lies; B.metanostictus, 1 $3 H.chinensis, 1 £; 
R limnocharis, 1%; R.guentheri. 1S. All of 
them are obtained from Geling, Yongtai, Fuzhou. 
1.2 Reagent 

Random oligonucleotide primer (Table 1), Taq 
DNA polymerase, 10 x PCR Buffer, 25 mmol/L Mg- 


Cl,. dNTP Mix (10 mmol each of dATP, dCTP, 
dGTP and dTTP), Rnase A, Proteinase K were from 
Shanghai Sangon Co. Other Reagents were made in 
China. 

Table 1 Random oligonucleotide primer sequences 


Primer Sequence Primer Sequence 
$301 CTGGGCACGA $311 GGAGCCTCAG 
$302 TICCGCCACC 5312 TCGCCAGCCA 
$303 TGGCGCAGTG $313 ACGGGAGCAA 
$304 CCGCTACCGA $314 ACAGGTGCTG 
$305 CCTTTCCCTC $315 CAGACAAGCC 
$306 ACGCCAGAGG $316 CTCTGTTCGG 
$307 GAGCGAGGCT $318 GACTAGGTGG 
5308 CAGGGGTGGA $319 TGGCAAGGCA 
3309 GGTCTGGTTG $320 CCCAGCTAGA 


5310 CCCTAGACTG 


1.3 Extracting and purification of genomic DNA 

Genomic DNA was extracted from the fresh tissue 
(liver, heart, muscle). After heing homogenated and 
lysis (0.1 mol/L Tris-HCl, al pH8.0, 5 mol/L Na- 
Cl, Smmol/L EDTA, 1% SDS, 100 ng/mL of pro- 
leinaseK, and 3h at 55°C), proteins and others im- 
purilies were removed by three sleps of phenol-chloro- 
form extraction, After an RNase (100 pg/mL, 1h at 
37°C) trealment was performed, the RNase was re- 
moved by phenol-chloroform extraction. DNA was then 
precipilated and resuspended in sterile distilled waler. 
Purity control of extracted DNA was assumed by mea- 
suring ralio of absorplion values al 260 and 280 nm. 
DNA concentration was evaluated using absorption val- 
ues at 260 nm. 
1.4 Amplification conditions 

Amplification reactions were performed in volumes 
of 20 uL containing 10 mmol/L Tris-Cl. pH9.0, 50 
mmol/L KCl, 2.5 mmol/L MgCl, 1% Nonidet P40, 
250 pmol/L each of dATP, dCTP, dGTP and dTTP, 
0.2 pmol/L primer, 50 ng of genomic DNA and 1.0u- 
nil of Taq DNA polymerase. Samples were subjected 
to PCR in a BIO-RAD Gene Cycler using 42 cycles of 
1 min at 94, 1 min al 37, 2 min al 72. 
1.5 Agarose gel eletropheresis 

Amplification products were analyzed on a 1.4% 
agarose gel with 0.5 x TBE buffer for 3 hours under the 
condition of 3 yoltage‘cm and room temperature. Then 


the eletropherograms were showed by staining with 
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ethidium bromide. photographed under ultraviolet 
lamp in BIO-RAD Gel Doc-1000. 
1.6 Data analysis 

Presence and absence of a given fragment that am- 
plified from genomic DNA were represented by “1” and 
“0” characters, respectively, Clear or weak and re- 


m 


productive bands were recorded as “I”. No or weak 
and none-reproductive bands were recorded as“0”. 
Thus a matrix of characters was obtained for phyloge- 
netic analysis as initial data. All the data were com- 
pared via RAPDistanse 1.04 software package by cal- 
culation of the NL coefficient (Nei & Li, 1979): 
NE = 2a/(b +c) 

Where a is the number of similar bands from two sam- 
ples, and b and c are the total number of bands from 
each sample. The NL coefficient allows a value of the 
similarity between two samples. This analysis provided 
data about the genomic distance between different sam- 
ples. The NL coefficient was used because it is recom- 
mended for routine computation of genetic similarities 
using RAPD data (Lamboy, 1994). The NE coeffi- 
cient gaye similarity values in the range 0 to 1 and were 
converted to genetic distances (D) as D = 1-NE. 
Cluster analysis was made for these 4 species and the 
systematic trees were built up respectively by UPGMA 
(Unweight Pair-group Method with Arithmetical Aver- 
age) method (Sneath & Sokal, 1973) and NJ (Neigh- 
bour-Joining ) method (Saitou & Nei, 1987) via 
MEGA 1.02 statistical software package . 


2 Result 


2.1 RAPD Result 

4 samples were screened using 1910-mer primers, 
16 except $307, S309 and 5310 with clear and repro- 
ductive finger-prints of interest generated. Fig. 1 shows 
profiles obtained with the $311 and $301 primers. 
2.2 Genetic distance and cluster analysis 

To show genetic affinity. the genetic distances a- 
mong these 4 species was calculated via RAPDistanse 
software package, as showed in Table 2. 

Fig.2 shows dendrograms that were build up re- 
spectively by UPGMA (Unweight Pair-group Method 
with Arithmetical Average ) method and NJ (Neigh- 
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bour-Joining) method via MEGA 1.02 statistical soft- 


ware package. 
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Fig.1 The electrophoresis patterns of RAPD of 4 species 

after random amplification by primers (a: prmer 

5311, b: primer 5301) 

A: B. melanostictus: B: H.chinensis; C: R, limnocharis; D: 
R. guentheri; M: molecular size marker (DNA ladder mit, MBI); 
N: negative control (H-0) 
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Fig.2 Dendrogram of 4 species 
a.UPGMA method; b.NJ method. 


3 Discussion 


Research on the population genetic and systematic 
evolution of Urodela in amphibian using RAPD method 
have been reported ( Wang et al., 1998) while the 
same kind of report on Anura lacked. In our research, 
checked by RAPD technique, we study the phyloge- 
netic relationship among B . melanostictus, H . chinen- 
sis, R.timnocharis and R . guentheri, which belong 
to the three differrent family of Anura. 


Totle numbers of the clear and reproductive of in- 
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Table 2 Genetic distances among 4 species 


Species B . melanostictus H. chinensis  R.hmnocharis R. guentheri 
B. melanostictus = 
H . chinensis 0.742331 一 
R . limnocharis 0. 798683 0.621918 = 
R. guentherr 0.911950 0.809524 0.577465 = 


terest bands generated by 16 primers, except 5307, 
$309 and $310 is 219. Polymorphism was showed in 
the banding pattern of these finger-prints. The numbers 
of bands generated by each primer fluctuate from 4 to 
25. Among all 219 bands obtained, 148 bands existed 
respectively in the profiles of each species, account for 
67.58%, appeared among the 4 species for only one 
time; 60 bands were shared by 2 species, account for 
27.39%; 10 bands were shared by 3 specie, account 
for 4.57%; only 1 band appeared in the profiles of all 
4 species, account for 0.47%. 

As the values of genetic distance in Table 2 
demonstrate that; The genetic distance between R. Zi- 
mnocharis and R . guentheri (0.577465), which are 
in the same family of Ranidae., is the nearest. The ge- 
netic distance between B . melanostictus and H . chinen- 
sts (0.742331) are closer than that between B . melan- 
ostictus and R.limnocharis (0.798883 ), 
B . melanostictus and R . guentheri (0.911950), be- 
tween H .chinensis and R . guentheri (0.809524) or 
that R . limnocharis 
(0.821918). The genetic distance between B . melan- 


between 


between H.chinensis and 


ostictus and R . guentheri is the furthest of all. The dif- 
ference of the genetic distance among the 4 species in- 
dicated that the relationship between Bufoidae and Hyli- 
dae is closer than that between Bufonidae and Ranidae 
at genomic DNA level. 

The cluster analysis also indicated that R . fimnoc- 
haris and R . guentheri ,which both belong to Ranidae, 
were Clustered together first; then B . meanosttictus and 
H . chinensis were clustered together; at last, all 4 
species were clustered together. The same result was ob- 
tained by the UPGMA and NJ methods. 

From above, it is confirmed that the relationship 
between Hylidae and Bufoidae is closer than that be- 
tween Bufonidae and Ranidae. Our result, at genome 
level, also provides new evidence of the conclusion 
about the systematic evolution of the three families of A- 
nura that drawed from the studies of morphology, chro- 
mosomal and mitochondrial DNA. Meanwhile, it also 
indicated that RAPD technique holds considerable 
promise for the population genetic, systematic evolution 


and classification of Anura. 
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WE: 参照 人 SRY 基因 HMG-box 保守 区 的 序列 ， 设 计 1 对 兼并 引物 ， 扩 增 了 乌 包 C Chinemys reevesii) 的 


Sos 基因 ， 并 对 扩 增 产物 进行 了 克 隆 和 测序 。 结 果 在 峻 雄 个 体 中 均 科 选 出 4 个 不 同 的 Sox 基因 ， 无 性 别 差 异性 ; 
其 DNA 序列 和 编码 的 氨基 酸 序列 与 人 相应 的 SON BER AY HAE HE BI 92%., 91% . 84%, 92% AT 100% ， 


93% , 98% 、98 允 ， 显 示 出 高 度 的 保守 性 。 


XRH: 乌龟 ; Sox 基因 ; 克隆 ; 测序 
中 图 分 类 号 : 9953, 9959.63 


SRY 基因 (sex-determining region of Y chromo- 
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用 RAPD 分 析 四 种 无 尾 类 的 系统 演化 关系 


林 玮 ARS KEK Hh H REE 高 建 民 
【福建 师范 大 学 生物 工程 学 院 福州 ”350007 hioforest @sina. com) 


HE. 利用 RAPD 技术 检测 了 分 属 无 尾 目 3 个 
科 CHASER. MRR AL. EPO KERERE (Buo 
melanostictus), P HAE ( Hyla chinensis), Y$ fti 4E 
(Rana limnocharis), WKH (R . guentheri) 的 系 
统 发 生 关 系 。 经 19 个 随机 引物 对 4 个 物种 基因 组 
DNA 进行 扩 增 ， 选 择 其 中 扩 增 谱 带 清晰 的 16 个 引 
物 进行 分 析 ， 计 算 不 同 科 间 及 同一 科 内 不 同 种 间 的 
HEEB. AREH: 16 个 引物 获得 的 RAPD 谱 
带 均 表 现 出 不 同 程度 的 多 态 性 ; 泽 陆 蛙 与 沼 水 蛙 间 


ASR i, MRES- pE EA 
SRK FR BER AE MS RE YP RE | rk ESL, 
及 中 国 雨 蛙 与 译 陆 蛙 、 沼 水 蛙 之 间 的 亲缘 关系 近 。 
从 基因 组 DNA 水 平 上 也 说 明 雨 蛙 科 与 昌 蜂 科 间 的 
亲缘 关系 更 近 ， 与 蛙 科 的 亲缘 关系 更 远 。 这 与 形态 
学 、 染 色 体 和 线粒体 DNA 多 态 性 研究 的 分 析 结 果 
一 致 ， 从 而 进一步 从 分 子 水 平 上 为 无 尾 目 这 3 科 的 
系统 演化 提供 了 新 的 证 据 。 


关键 词 : 随机 扩 增 多 态 DNA: 无 尾 目 ; 系统 演化 ; 基因 组 DNA 
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